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Abstract 

This paper reports measurements using the technique of On Line Nuclear Orientation (OLNO) 
which reexamine the gamma decay of isomeric ^^OHP" and specificahy the 501 keV 8~ - G"*" tran- 
sition. The irregular admixture of E2 to M2/E3 multipolarity in this transition, deduced from 
the forward-backward asymmetry of its angular distribution, has for decades stood as the prime 
evidence for parity mixing in nuclear states. The experiment, based on ion implantation of the 
newly developed mass-separated i^Ojjpi ijg^m at ISOLDE, CERN into an iron foil maintained at 
millikelvin temperatures, produces higher degrees of polarization than were achieved in previous 
studies of this system. The value found for the E2/M2 mixing ratio, e = -0.0324(16) (17), is in 
close agreement with the previous published average value e = - 0.030(2), in full confirmation of 
the presence of the irregular E2 admixture in the 501 keV transition. The temperature dependence 
of the forward-backward asymmetry has been measured over a more extended range of nuclear 
polarization than previously possible, giving further evidence for parity mixing of the 8~ and 8"*" 
levels and the deduced E2/M2 mixing ratio. 

PACS numbers: 21.10.HW,23.20.En,23.40.Bw,27.70.+q,29.30.Lw,150<A<189 
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I. INTRODUCTION 



Parity, reflection symmetry in the origin of a co-ordinate system, is one of the fundamen- 
tal symmetries of physics. Establishing whether parity may be taken as a conserved quantity 
under a system of forces, or alternatively determining the conditions under which, and the 
degree to which, parity is not conserved, form basic constraints upon physical theories. The 
discovery of parity non-conservation (PNC) in the weak interaction was one of the most 
important discoveries of modern physics. However the extent to which parity is to be con- 
sidered a conserved quantity in nuclear phenomena remains a challenge to both experiment 
and theory. Parity mixing in bound nuclear systems is understood as a consequence of weak 
(parity violating) interaction terms in the nuclear Hamiltonian and precise calculations of 
this phenomenon are not yet available. 

Of the many experiments aimed at detecting parity non-conservation in nuclear states, 
one, the measurement of an irregular E2/M2 mixing in the 8~ - 6"^, highly K-forbidden, 
501 keV gamma decay of the 5.47 h isomer of ^®°Hf, stands out. It was first observed in 

riQ 

gamma ray circular polarization experiments |2j . Using the technique of low temperature 
nuclear orientation (LTNO), the reported result is an apparently well established effect, a 
32/M2) = -0.030(2), of a magnitude more than fifteen times the experimental 
error U, y, gI. Other LTNO measurements are listed in [7]. Two other statistically 
significant PNC results, of much smaller experimental effects, on states in ^'^^Lu [sl and ^'^^Ta 
[9|], have been reported using circular polarization technique. On no other bound nuclear 
systems, including the recent work on ^^^j^m ^^iq latest published reports claim to 

:ind any effect deviating by more than two standard deviations from parity conservation 
3, lH. The result on the i^ogpi jgomeric decay, based on experimental evidence of both 
angular distribution studies from nuclei polarized at millikelvin temperatures and circular 
polarization studies (see p] and references therein), dating from the 1970's, stands today as 
the strongest evidence of the level to which nuclear states cannot be taken as eigenstates of 
parity. 

It is important to examine the validity of this significant result. This paper describes 
a remeasurement of the evidence for parity non-conservation effects in the ^^''Hf isomeric 
decay, using the techniques of on-line nuclear orientation available today. Developments of 
technique over 30 years have given access to different source making methods, leading to 
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higher degrees of polarization and the abihty to make continuous observations over a period 
of days rather than successive measurements on a series of decaying samples. 

The paper starts with description of previous work, a brief introduction to the necessary 
formalism and the justification of a new study in Sec. [TTl This is followed (Sec. Illip by 
the detailed description of the new experiment and the analysis of the gamma ray spectral 
data. Comparison with theoretical calculation of the observed effects requires discussion 
of aspects of the hyperfine interaction and of the angular distribution coefficients, given in 
Sec. HVl Sec. |V]includes discussion of the quality of implantation and aspects of thermometry, 
leading to evaluation of the parity violating E2/M2 mixing ratio in the 501 keV transition. 
A brief discussion of the final results and the need for, but difficulty of, its interpretation in 
nuclear theory, is given in Sec. IVII 

II. FORMALISM, PREVIOUS EXPERIMENTS AND TECHNICAL DEVELOP- 
MENTS. 



Full descriptions of the formalism of the nuclear orientation experimental method and 
of the previous measurements on i^oj^pi ^^^.^ given in Refs. Here only the 

essential details relevant to the present work and its analysis are presented. Fig [1] shows a 
partial level scheme of ^^°Hf with the 5.47 h, V= 8~ , K=8, isomer decaying to the 8"^ and 
G"*" levels of the K=0 ground state band. In a nuclear orientation experiment the radioactive 
sample is polarized by cooling to millikelvin temperatures in an environment in which the 
nuclei experience a magnetic field sufficient to produce a high and controllable degree of 
nuclear polarization. The gamma radiation from such a polarized source exhibits a strongly 
anisotropic angular distribution. If the nuclear states between which the gamma decay 
occurs are eigenstates of parity, the theoretical distribution contains only even terms in 
cosine theta, measured from the axis of polarization, and is given by the expression loj] 

W{e,T) = 1 + E fBx{T)UxAxQxPx{cose) (1) 

A even 

where Ba(T) are the temperature dependent orientation parameters, describing the polar- 
ization of the parent nuclei, Ua are angular momentum coupling constants describing the 
change in the polarization when the decay passes through a sequence of unobserved transi- 
tions between states before the detected gamma emitting state is reached, A^ are angular 
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FIG. 1: Decay of the 5.47 h i«ujjpn 

momentum coupling constants describing the multipolarity of the observed transition, Qx 
are sohd angle correction factors for finite detector size and Pa are associated Legendre 
polynomials. The factor /, the fraction of nuclei in sites which experience the hyperfine 
field, is particular to implanted sources and is fully discussed in Sec. IV A[ It is usual to 
place detectors at 0°, 90° and 180° to the polarization axis and to describe the orientation 
in terms of the anisotropy a(T) given by 

W{0°,T) 



a{T) 



1. 



(2) 



1^(90°, T) 

The normalized intensities at 0° and 180° are the same under the assumption of parity 
conservation. 

If either of the nuclear states between which the gamma decay occurs are not eigenstates 
of parity, odd A terms in the distribution, of the same form as Eg. ^ are introduced. These 
terms are seen most clearly in the asymmetry ^(T), defined in 

Vr(0°,T) - Vr(180°,T) 



I dl^ 
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V(0°,T) + 1^(180°, T) 



2 E fBx{T)UxAxQxPx{cosd) 

A odd 

1+ E fBx{T)UxAxQxPx{cos( 

A even 



(3) 
(4) 



The previous work sought evidence for parity non-conserving admixtures in both the 57 keV 
8^ - 8+ transition and the 501 keV 8^ - 6+ transition by measuring the asymmetry between 
the gamma intensity measured at zero degrees and 180 degrees to the axis of polarization of 
samples of ^^ogpn (HfxZri_x)Fe2 alloys cooled to millikelvin temperatures 

mm 

. For 

the 57 keV transition the reported gamma ray asymmetry was small, consistent with zero 
to one part in a thousand. However the 501 keV transition consistently showed asymmetry 
of close to 1.5% at a temperature of about 20 mK, more than ten times the experimental 
error. 

Two, related, major technical developments over the thirty years since the last reported 
experiments on i^ogpi ^^gQ^y make an experimental re-examination of the case timely. The 
first of these is the development of the on-line nuclear orientation technique in which a beam 



of radioactive ions is implanted into a ferromagnetic metal host foil cooled to mi 



temperatures in a ^He/^He dilution refrigerator with side access for the beam 12(]. This 



likelvin 



method allows continuous measurement at a controllable range of temperature down to below 
15 mK with a uniform source strength. In contrast, the previous experiments involved a 
series of samples cooled by contact with a demagnetized paramagnetic salt, which thereafter 
steadily warmed giving a slowly changing degree of nuclear polarization and each successive 
sample decaying with the 5.47 h half-life giving rise to variable dead-time and electronic 
pile-up correction. The second is the development of a means to produce a beam of ionized 
Hf, a highly refractory element, from an isotope separator. 

As well as gaining access to better controlled temperatures and a steadier source strength, 
the use of ion implantation allows use of pure iron, rather than the (HfxZri_x)Fe2 alloy, in 
which to polarize the nuclei. The hyperfine field at Hf in iron is close to 70 T (see below), 
more than 3 times stronger that the ~ 20 T present in the alloy. This means that the hafnium 
nuclei can be more fully, indeed almost completely, polarized at attainable temperatures. 

In view of these basic improvements of available technique a new experiment on the 
anisotropic emission of gamma radiation from oriented i^oj^pi j^g^g Iqqqt^ carried out. 
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III. THE EXPERIMENT 



A. Outline 

The experiment was performed at the ISOLDE isotope separator facihty, CERN, using the 
NICOLE on-line nuclear orientation dilution refrigerator system [l^. The 1.4 GeV protons 
were incident upon a mixed Ta/W metal foil target and the hafnium atoms produced were 
transfered to the ion source using fluorine, added to the plasma support gas as CF4. The 
most intense emerging hafnium fluoride molecular ion beam, used for the experiment, was 



HfFj^ [13|. The Hf orientation sample was prepared by accelerating the ions to 60 keV and, 
after mass separation, impinging them on the surface of a pure (99.99%) iron foil soldered to 
the copper cold finger of the dilution refrigerator, which is perpendicular to the separator ion 
beam (see Fig. [2]). At impact the molecular ions disintegrated and the hafnium nuclei were 
implanted into the iron lattice, producing a source entirely free of contaminant activity. The 
F nuclei also enter the iron but come to rest in a region well removed from the Hf nuclei. 
The total implantation dose was approximately 4 x 10^^ ions into a region of ~3 mm in 
diameter giving a maximum local concentration of Hf in the iron foil below 10^^ atomic 
percent. The quality of the implantation, as indicated by the fraction of hafnium nuclei 
stopping in substitutional lattice sites in the iron, is discussed below. A second iron foil, 
containing diffused ^''Co activity for which all orientation properties are known, was soldered 



to the back of the cold finger to act as a nuclear orientation thermometer iJ]. During the 
experiment the iron foil sample was magnetized to saturation using a Helmholtz pair of 
polarizing coils which produced a magnetic field of 0.5 T applied in the plane of the foil 
and at right angles to the beam axis. The direction of the field, which determines the sense 
of polarization of the radioactive nuclei, could be reversed by changing the direction of the 
current through the coils, a procedure which took about 10 minutes. 

Gamma radiation from both activities was detected in three large intrinsic germanium 
detectors, two placed on the field axis on opposite sides of the cryostat, at 0° and 180° to the 
axis of polarization (depending upon the field direction) and one at 90° to the axis of polariza- 
tion below the cryostat. The spectra were accumulated in files of 300 s duration throughout 
the experiment, pulses from a regular pulser being introduced to the pre- amplifiers to allow 
correction for pile-up in detectors and dead time in the electronics. 
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FIG. 2: Schematic experimental set-up. 



The gamma ray spectrum (see Fig. [3]) contained six strong, fully resolved, transitions, at 
501 keV, 444 keV, 332 keV and 215 keV from is^Hr (the 57 keV and 93 keV transitions 
were strongly absorbed in cold finger and cryostat), and at 137 keV and 122 keV from 
the ^^Co thermometer. To obtain counts for each transition detected in each 300 s file, 
simple windows were set over the peaks, taking care to set them wide enough so that any 
small changes in gain of the system did not result in loss of counts at any time during 
the experiment. Background to each peak was found by setting additional windows on 
the spectra above and below the peaks and making a linear interpolation to determine the 
background count in each peak window. 

The sample was cooled to achieve high degrees of nuclear polarization and the polarizing 
field reversed at intervals when the refrigerator was operating steadily and the implanted 
ion beam was stable. At each reversal the change in gamma intensity recorded in each of 
the two axial detectors was measured in all transitions. For transitions in which parity is 
conserved there should be no change in intensity since any observed change is proportional 
to the parity non-conserving transition amplitude. The principal objective of the experiment 
was to observe the temperature dependence of the parity non-conserving amplitude reported 
in the 501 keV transition over a wider range of temperature, and hence to higher degrees of 
nuclear polarization, than had been accessible to the previous studies (see Fig. Hj). 

After the beam was first introduced to the refrigerator, the i^ogpn activity was allowed to 
accumulate for about 5 hours with the temperature close to 1 K and the iron foil unpolarized. 
The polarizing field was then applied and a sequence of 'warm' reference spectra were taken 
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FIG. 3: Gamma-ray spectrum showing isOjjpn transitions at 501, 444, 332 and 215 keV and ^'''Co 
transitions at 137 and 122 keV, with no Pb absorber present. 

with negligibly small nuclear orientation. As the source strength approached its asymptotic 
value, the sample was cooled to about 13.5 mK, the lowest accessible temperature with the 
beam present. Just after this, interruption of the beam from ISOLDE reduced the heat 
input, allowing further cooling to the refrigerator base temperature, 7.6 mK. This complete 
sequence is referred to as the initial cool-down. Shortly after the beam had returned, the 
first of a series of eight reversals was carried out. For each reversal the temperature was 
measured using nuclear orientation thermometry as described later. The first six reversals 
were done with the implantation beam incident on the sample. They were accordingly at 
temperatures which reflected balance between beam and radiation heating to the cold flnger, 
plus heating from absorption of the gamma activity in the source, and the cooling power 
of the refrigerator. The temperatures ranged between 15.6 and 25 mK, increasing slowly 
with time as the beam intensity and source strength grew. Later the mixing chamber of the 
refrigerator was heated to give a temperature close to 60 mK, still with implanted beam, 
and a sixth reversal took place at this relatively high temperature. Finally at the end of 
the experiment, when implantation was stopped, the sample cooled to the refrigerator base 
temperature of 7.6 mK during the decay of the activity and the seventh, lowest temperature, 
reversal was made. An eighth, 'dummy', reversal was made after the sample was warmed to 
about 1 K when no orientation was present. 
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B. A first survey of the results. 

For the 501 keV transition to show appreciable asymmetry, assuming that the other tran- 
sitions in the isoj^pi ^jg^a^y qj-q symmetric, a simple way to show the effect is to consider the 
ratio of the peak counts in the (asymmetric) 501 transition to those in another (symmet- 
ric) transition, for example that at 444 keV, or combination of symmetric transitions. As 
is shown in Appendix A (for the case in which the sum of the 444 keV and 332 keV counts 
is taken as the symmetric 'norm') the change in such a ratio, measured in a single detector 
on the axis of orientation, when the direction of nuclear polarization is reversed, is directly 
proportional to the asymmetry of the 501 keV transition. 

Fig. m (top panel) shows the measured double ratio 

^ iV(501,L)/iV(444,L) 

N{501,R)/N{AAA,R) ^ ' 

for each file during the experiment. Here L and R refer to the two detectors on the axis of 
nuclear polarization. When the applied field, and hence the direction of nuclear polarization, 
is reversed, the asymmetry, and with it any difference in this ratio from unity, must change 
sign. Fig. m (bottom panel) shows the same ratio, but for the 332 keV transition to the 444 
keV transition. The file numbers at which field reversals were made are indicated by the 
sequence of upward and downward arrows. In both panels the scatter of the data indicates 
the statistical error. It is immediately apparent that the data in the upper panel show 
variation about unity which reverses at every field reversal except the eighth. This last 
reversal took place above 1 K, with zero nuclear polarization, so a null effect is expected. 
By contrast, the data in the lower panel exhibit no such field direction dependence, as is 
expected for transitions between good eigenstates of parity. It is however also clear that the 
data in the upper panel show, in addition, some variation from unity when the temperature 
of the dilution refrigerator, i.e. the degree of nuclear polarization and hence the count 
rate in the different detectors, is changing rapidly. Cooling and heating took place over 
relatively narrow ranges of file numbers whose centres are indicated in the top of the figure. 
The apparent asymmetry suggested by these additional changes is due to different response 
of the detectors to strongly changing count rate. Such spurious effects have been fully 
eliminated in the more comprehensive analysis presented in the following sections. 
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C. Detailed description of the experiment and spectral analysis of the field rever- 
sals. 

In this section the procedure is described in sequence and in considerable detail, allowing 
full discussion of the several stages of the experiment. 

1. Selection of data files for analysis. 

The sample temperature at any time is determined by a balance between a combination 
of heating and cooling processes. During the early part of the experiment many short 
interruptions of the beam took place, during which the sample temperature fell rapidly over 
even just a few 300 s data files, recovering to its pre-interruption value within a few files 
of the return of the beam. All files involved with such variations of temperature were set 
aside, as were others during which the dilution refrigerator was being filled with cryogens, 
which can also somewhat perturb its operation. The files used in the analysis were free of 
any detectable outside perturbations. During each magnetic field reversal the current in the 
polarizing coils was first reduced over a period of about 300 s then the leads to the power 
supply were manually interchanged and the current increased over a similar period. This 
procedure resulted in small changes in the refrigerator temperature as detected in all gamma 
transitions, caused by eddy current heating in the cold finger holding the sample and the 
mixing chamber. The initial temperature was recovered in a further three or four files. Data 
files taken during these periods of heating and recovery were also discarded in the analysis 
presented below. 

2. The measured quantities. 

Since the source strength during implantation is a variable quantity, all anisotropy and 
asymmetry measurements arc in the form of ratios of counts in detectors at different angles 
to the polarizing field. Furthermore, to eliminate effects of variable dead time and pile up 
during different counting periods, each raw peak count is divided by the pulser peak count 
for the same file and detector, referred to as the pulser normalized count. The quantity 
Wexp(^) is the ratio of the pulser normalized count in a specific gamma peak in a detector 
at angle 9 to the polarization axis, to the pulser normalized counts in that peak from an 
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unpolarized sample, the 'warm' counts. The measured value of anisotropy a{T) is given by 
the ratio [Wexp(0°)/Wexp(90°) - 1] from which temperatures are deduced and the measured 
asymmetry ^(T) is given by 2[Wexp(0°) - Wexp(180°)]/[Wexp(0°) + Wexp(180°)], sensitive to 
the degree of parity non-conservation. Both have been defined above. 

3. Introduction of absorbers. 

The initial cool-down, combined with the later cool-down after the implanted beam was 
finally stopped, provided data allowing comparison of the anisotropies measured on the 
isojjpi ^''Co transitions over a wide temperature range. Since temperatures can be 
deduced from the observations on ^^Co this allowed calibration of the i^OHf™ transitions as 
secondary thermometers so that counts from the ^^Co transitions were no longer required 
for temperature measurement. Lead absorbers of equal thickness, giving rise to attenuation 
factors (measured to be equal to ±1%) for a given gamma energy in all detectors, were 
secured over each detector face to reduce the total counting rates. The close equality of the 
attenuation factors allow use of the pulser normalized 'warm' count ratios in analysis of the 
ratios measured with absorbers present. 

4- Field reversals 1-6. 

Following the initial cool-down, once the beam returned and the temperature had reached 
equilibrium, the first field reversal was carried out. Further reversals were made at intervals 
of about 7 hours, the time between them being determined largely by the requirement 
of steady beam and good temperature stability. The results of all field reversals and the 
conditions under which they were performed are summarised in Table [H During this sequence 
the temperatures deduced from the i^opjpi transitions rose slowly from 15.6 mK to close 
to 25 mK as the sample strength grew, reflecting a slow increase in yield from the ion 
source. Reversals 2-4 were carried out in the same experimental set-up. For reversal 5 
the pulser, which produced some distortion of the gamma transition peak shapes through a 
small undershoot in the amplifier base line, was removed to check that its presence did not 
affect the measurements. Lack of pulser normalization for this reversal (and for reversal 7) 
meant that rather than being able to evaluate the asymmetry for each transition separately. 
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only the difference of asymmetry between a pair of transitions can be extracted. This is 
discussed fully in Sec. IV CI and in Appendix A. 

The pulser was reintroduced before reversal 6, for which the refrigerator temperature was 
raised to close to 60 mK by supplying heat to the mixing chamber. For this reversal ^^Co 
thermometry was needed, thus the lead absorbers were removed. 

5. After-beam cool-down, low temperature field reversal 7 and final warm-up. 

The isojjpn ]3ga^]2i was stopped shortly after reversal 6 and the dilution refrigerator cooled 
fairly quickly to its base temperature as the sample slowly decayed. The absence of lead 
absorbers during the cooling allowed use of the ^^Co thermometer. Reversal 7 was carried 
out when the temperature had reached its equilibrium value of 7.6 mK and was made with 
lead absorbers present but without use of the pulsers. When this reversal was complete 
the refrigerator was warmed to above 1 K and a 'dummy' reversal was performed whilst 
additional 'warm' reference files were recorded for approximately 24 hours. The ratios of 
pulser normalized gamma peak counts from the initial and final 'warm' files in all detectors 
were unity to within 0.6%, demonstrating that there had been no significant drift. It should 
be made clear that this 0.6% variation over the course of the experiment does not affect the 
asymmetries derived from the data. This follows since upon each reversal the asymmetry is 
the immediate change in the ratio of the normalized count rates in the 501 keV transition 
in the 0° and 180° detectors, so normalization to warm counts cancels. In interpreting 
the asymmetry using Eqn. H] the numerator is directly given by the reversal data. The 
denominator (the sum of the anisotropics) does depend on the 'warm' counts. However, in 
the even terms, a small change, such as the 0.6% referred to, simply adds to the anisotropy 
in one direction and subtracts from the other, leaving their sum unchanged. 
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IV. PARAMETERS USED IN CALCULATION OF THE ANGULAR DISTRIBU- 
TIONS. 



A. The magnetic moment of ^^°HP^. 

The 8~ isomeric state of ^sogpn gxpected to be a rather pure two-quasi-particle state, 
very similar to the 8~ isomers found in other Hf isotopes. Only one of these, I'^'^f^fm i^^g 
a measured magnetic dipole moment, //(^"^^Hf™, 8~) = 7.86(5) n.m. This moment, 

however, is in excellent agreement with predictions for a pure two-proton excitation com- 
prising g7/2+ and hg/2- protons. Taking the g-factors for these two protons from neighboring 
isotopes, simple combination gives prediction for the moment of the 8~ isomer between 7.7 
and 7.9 n.m. Although there are many measurements of the g7/2+ single quasi-particle state 
in the vicinity, and they show only a small variation with neutron number, there is only one 
accurate measurement of the h9/2- state (in ^^iTa) Q, thus any changes in the magnetic 
moment of the 8~ Hf isomers with addition of neutrons cannot be argued from experimental 
evidence. However they are likely to be small. 

Low precision measurements for the magnetic moment in i^ogpi i^^yQ been reported 
based on the hyperfine splitting extracted from LTNO in the Hf alloy samples used in that 
work (see below), giving results of close to 9 n.m. with errors of 1 n.m. These appear high 
compared with estimates based on the nuclear structure of the isomer. 

B. The hyperfine field acting at isohF" in iron. 

The magnetic hyperfine interaction at a substitutional site in cubic iron is a well defined 
quantity which involves the product of the magnetic dipole moment of the isotope concerned 
and the hyperfine field acting at that site. Until recently there were no accurate measure- 
ments of the magnetic dipole moment of any radioactive isotope_of hafnium [l^ and thus 



only inaccurate values of the hyperfine field could be obtained [17 [. However, recently the 
magnetic moment of ^''^Hf was measured to be - 0.677(9) n.m. 18| and nuclear magnetic 
resonance of ^^^Hf oriented in iron was reported at a frequency of 139.0(1) MHz These 
results yield the hyperfine field for Hf in iron as - 67.5(9) T. 

Estimate of the interaction strength for is^Hf™ in iron can be obtained from published 
experimental results only through taking the hyperfine interaction expressed as a tempera- 
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FIG. 5: Calculated Ba coefficients vs 1/T for hyperfine interaction 25.9 mK, each normalized to 
its saturation value at T = 0. The vertical dotted lines show the range of Ba accessed by previous 
measurements 3, 0, 0, 0]; full lines indicate the range covered by present work. 



ture, Tint= /iBhf/Ik = - 8.2(2) mK determined for i^oj^pn ^j^^ compound (Hfo.iZro.9)Fe2 
and the ratio 3.9(4) of the strengths determined from Mossbauer studies on the 93.3 keV 
state in iron, i.e. Tint = - 7.6(1) mK, and in the compound (Hfo.iZro.9)Fe2, Tint = - 1.93(17) 
mK {20I]. These values give the interaction in iron as Tint = - 32.4(30) mK. 

However, taken with the new, accurate, result for the hyperfine field, such an interaction 
would predict a magnetic moment of 10.4(10) n.m. for ^'^''Hf™, far larger than the measured 
moment in I'i'^fjfm ^^d above any reasonable increase due to additional neutrons in ^^°Hf™. 
This is considered to be too high and therefore an estimated interaction. Tint = - 25.9(13) 
mK, based on the measured hyperfine field and a magnetic moment of 8.4(5) n.m. consistent 
with, but somewhat larger than, the measured value in i72jjpi ^^^^ been used to calculate 
the orientation parameters Ba for isopjpn iron. The results for A = 1 - 6 are plotted in 
Fig. [5l as a function of inverse sample temperature 1/T. In the figure, each Ba is normalized 
to its saturation value. The range of the orientation parameters at which field reversals were 
carried out in this work is indicated in the figure and compared with the range accessed 
by the experiments of Refs. [3, 0, 0, O] showing the increase in the degree of polarization 
and the close approach to full saturation achieved in the present experiment. The very high 
degree of polarization achieved has the important consequence that, with uncertainty in the 
hyperfine interaction estimated to be 6%, the uncertainty in the extracted irregular parity 
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admixture e is less than 1%, small compared with the statistical errors. 



C. The nuclear spin-lattice relaxation time of ^^'^Hf'^ and Co in iron. 

To assure full thermal equilibrium between the iron lattice, taken here to be at constant 
temperature, and the implanted hafnium nuclei prior to their decay it is necessary that the 
hafnium nuclear spin-lattice relaxation time Ti be short compared to the nuclear lifetime. 
For nuclei in metals, relaxation is via the conduction electrons and can be described in 
terms of the Korringa constant = TiT, where T is the absolute temperature. For a given 
system in iron the value of Ck can be estimated from the empirical relation CkT 

int — 1.4 



X 10~^ sK^ 2l|] where Tint is the nuclear interaction strength. Tint = - 25.9 mK for ^soj^pn^ 
This yields Ck ~ 0.21 sK and an estimated relaxation time of 2.0 s at 100 mK. At lower 
temperatures the simple inverse relation to lattice temperature breaks down and Ti reaches 



a maximum value approximately given by 3.3 Ck/(I + |)Tint |22|, estimated for isoHf™ to 
be ~ 3 s, far shorter than the lifetime of 5.47 h, thus clearly satisfying the requirement for 
thermal equilibrium. 

As described below, knowledge of the fraction, /, of hafnium nuclei in substitutional 
lattice sites of the iron foil, is necessary both for analysis of the measured asymmetry of the 
501 keV transition and also to calibrate the hafnium gamma transition anisotropies for use 
as thermometers for part of the experiment. The determination of the fraction is done by 
comparing the anisotropies of the i^^Hf™ and ^^Co gamma transitions over as wide a range 
of temperatures as possible. When the lattice temperature is changing, for example during 
cool-down, for such comparison to be valid it is necessary that for both the hafnium and 
cobalt isotopes the spin-lattice relaxation times be substantially shorter than the time of 
measurement for one file. For is^Hf™ this is already established. For ^''Co Tint is 14.2 mK 

n 

and the Korringa constant is Ck = 0.4 sK [23] so that the relaxation time varies between 4s 
at 100 mK and a maximum of 23 s well below 14 mK, again both much shorter than the file 
measurement time of 300 s. Thus at all temperatures in this experiment the two isotopes 
can be taken as being in good thermal equilibrium with the iron lattice and with each other. 
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TABLE I: Analysis of both i^OhP' and Co gamma ray data to give asymmetry values for all field 



reversals. 



Reversal 




Pb 


Field 


Pulser 


T 


Asymmetry ^(T) (%) 


number 


Beam 




Direction 




(mK) 


122 keV 137 keV 


215 keV 


332 keV 


444 keV 


501 keV 


1 


ON 


OFF 


L - R 


ON 


15.6(2) 


+0.08(11) -0.26(47) 


+0.12(15) 


-0.07(14) 4 


-0.21(15) 


-1.09(32) 


2 


ON 


ON 


R - L 


ON 


19.2(13) 


lead absorber present 


-0.27(19) 


-0.09(10) 4 


-0.01(10) 


-1.24(19) 


3 


ON 


ON 


L - R 


ON 


20.8(16) 


lead absorber present 


-0.01(19) 


+0.01(10) +0.11(10) 


-1.12(19) 


4 


ON 


ON 


R - L 


ON 


22.2(14) 


lead absorber present 


+0.15(18) 


-0.12(10) 


-0.23(10) 


-1.41(23) 


5 


ON 


ON 


L - R 


OFF 


25.0(13) 


evaluated ratio 501/(444-1 


-332) see text 




-1.19(13) 


6 


ON 


OFF 


R - L 


ON 


57(7) 


-0.41(16) -1.89(73) 


-0.12(11) 


-0.06(11) 


-0.15(11) 


-1.48(26) 


7 


OFF 


ON 


L - R 


OFF 


7.6(1) 


evaluated ratio 501/444+332) see text 




-0.93(13) 


8 


OFF 


ON 


R - L 


ON 


>1000 


-0.11(6) -0.33(25) 


+0.14(8) 


-0.10(8) 


-0.13(8) 


+0.11(18) 


Weighted average asymmetry for 
















transitions other than the 501 keV 






-0.08(9) -0.74(74) 


-0.04(7) 


-0.07(5) 


-0.03(8) 




for cold reversals 1 


- 7. 



















D. The gamma transition angular distribution coefficients. 

Another measured parameter required in the angular distribution calculation is the 
E3/M2 multipole mixing ratio in the 501 keV transition. The value 5(E3/M2) = + 5.3(3) 
given in Ref. ^] has been taken. The uncertainty in this parameter produces an error in the 
final result for the irregular E2/M2 mixing ratio e again much smaller than the statistical 



errors. The 57 keV 8^ - 8+ transition was taken as pure electric dipole js]. All other tran- 
sitions are of pure electric quadrupole multipolarity. The intensities of the 501 keV and the 
(unobserved) 57 keV - 444 keV decay paths feeding the 641 keV 6"^ level were taken as 14% 



and 86% of the total respectively |24|. 



In Table HTIthe calculated values of the even term angular distribution parameters UaAaQa 
for each analysed transition in the decay of i^oj^pn given. The parity non-conserving 
terms for the 501 keV transition are given to first order in the irregular mixing ratio e by ^ 

= T^[^a(2288) + 5Fa(2388)]. (6) 

For e = - 0.030 and 5 = + 5.3 the values are: UiAiQi = - 0.0045, U3A3Q3 = 0.0016 and 
U5A5Q5 = 0.0019. Here the value of epsilon from [6] is taken; the variation of the odd A^ 
with epsilon is given by Eqn. [6] and illustrated in Fig. [HI 
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V. FURTHER ANALYSIS. 



A. The fraction in good sites: calibration of the i»OHf™ thermometer. 

One non-ideal feature of the use of ion implantation for sample preparation is the fact 
that a fraction of the implanted hafnium nuclei come to rest at sites which do not experience 
the full, substitutional site, hyperfine interaction. Such sites include not only irregular sites 
in the iron matrix, but also nuclei which undergo strong interactions in the surface layers of 
the target and come to rest in the thin oxide layer which is always present on iron foils. The 
low-temperature nuclear orientation technique does not have the ability to explore details 
of the site distribution, but investigations have shown that frequently a valid description of 
the system is to consider a two-site model, with fraction / in the substitutional site and the 
remainder in a zero field site thereby remaining unoriented at all temperatures. The model 
can be shown to be valid in a particular case if the fraction / extracted using it is found to 
be constant, independent of the sample temperature. 

The substitutional site hyperfine interaction is known to about 6%, as discussed above, 
providing an extremely high degree of polarization of this fraction at the lowest temperatures 
reached in this work. The angular distribution Wcai of the i^oj^pn g^j^ima transitions from 
nuclei in such sites can be calculated as a function of temperature from their known multipole 
character. A value of the fraction / can then be extracted from the data by comparison of 
the measured anisotropy Woxp of each of the 215, 332, 444 and 501 keV transitions, provided 
the temperature is known, using the relation [Wexp(6') - 1] = f [Wcai(^^) - 1]- 

The sample was cooled from above 50 mK to the base temperature of 7.6 mK twice during 
the experiment, the initial and final cool-downs, both times with temperatures determined 
from the observed anisotropics of the pure electric quadrupole 137 keV transition of ^^Co. 
The isoHf™ data on all gamma transitions from all files taken during both cool-downs and 
at base temperature have been analysed for the fraction /. The results are shown in Fig. El 
It is seen that a value of / close to 80% is obtained from each transition, over the full 
temperature range of the experiment, with only a slight upward drift (1-2%), justifying use 
of the simple two-site distribution model. The value adopted for further analysis of the data 
is / = 0.805(10), as indicated in Fig. [61 

With the fraction / determined, measured anisotropics of the i^oj^pn gg^j2ima transitions 
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FIG. 6: Fraction in good sites, /, evaluated from data on i^Ojjpn ^jg^ay transitions as described in 
the text. The full lines give the range of the determined value /=0. 805(10). 



can be used as thermometers. This was done in analyzing data of reversals 2,3,4, and 6, 
when the presence of lead absorbers removed the ^''Co transitions from the spectra. 



B. Temperatures of reversals 5 and 7. 

These two reversals were carried out without the pulser peak present in the gamma 
spectra. This means that whilst relative anisotropies and asymmetries of different gamma 
transitions can be evaluated, absolute, source strength corrected, values for individual tran- 
sitions are not available and thus temperature was not directly measured. However indirect 
arguments can be made to establish the temperature for these reversals. 

isojjpi transitions with pulser present were available for thermometry until shortly before 
reversal 5 and after it the pulser was reconnected and lead absorbers removed so that both 
isojjpi g^^^ Co thermometry were available. At this stage of the experiment the implanted 
beam was steady with time and the sample activity had reached its asymptotic value, thus 
there was no reason to expect the sample temperature to vary. Fig. [7] (upper panel) shows the 
ratio of the number of counts in the 444 keV transition peak recorded in one axial detector 
to the number recorded in the 90° detector during the period containing the reversal. The 
figure shows constancy of this ratio, that is constancy of the sample temperature, to within 



20 




540 



680 700 720 740 

Data acquisition file number 



FIG. 7: Upper panel: Ratio of counts in the 444 keV gamma transition peak in one axial (0°, 180°) 
detector and the 90° detector for the time period spanning reversal 5. The generally constant level 
of this ratio shows constancy of temperature and determines the temperature of reversal 5. The 
variations seen between files ~ 430 and 470 are caused by interruption in the hafnium beam, whilst 
those around file 485 are due to the field reversal itself. The ranges of files used for pre-reversal 
and post-reversal averages are indicated by vertical arrows. Lower panel: As in the upper panel 
but for the period spanning reversal 7, showing the data used to determine the temperature of 
reversal 7 (see text). The step at file 670 was produced by dead-time changes when the pulser was 
removed. The statistical errors are smaller than the data point symbols. 



about one millikelvin over this period. ^''Co thermometry shortly after the field reversal 
(files 540-548) gave the temperature as 25.0(13) mK, which is assigned to reversal 5. 

Reversal 7 was performed about 6 hours after the beam was stopped. The dilution 
refrigerator cooled sharply, reaching base temperature at least 3 hours before the reversal. 
The 444 keV un- normalized 0°/90° ratio as above is shown in Fig. [7| (lower panel), again 
indicating constant temperature after the cool-down over the period of the reversal. The 
^"^Co temperature measured at the end of the cool-down was 7.6(1) mK, which is taken as 
the temperature for reversal 7. 

The 'dummy' reversal 8 was performed after warm-up of the source to above 1 K and 
analysed as the other reversals; as required the average asymmetry was zero within statistical 
errors for all transitions (see Table H]) . 
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TABLE II: Angular distribution coefficients used in calculation of the even polynomial terms in 
the anisotropy of transitions from ^^'^Hf™. 



E^ [keV] 


ir 


TTT 


Multipolarity 


U2A2Q2 


U4A4Q4 


UeAeQe 


501* 


8- 


6+ 


M2/E3 


-0.407(11) 


-F0.106(9) 


-0.222(8) 


444 


8+ 


6+ 


E2 


-0.338(1) 


-0.114(1) 


0.000 


332 


6+ 


4+ 


E2 


-0.332(1) 


-0.107(1) 


0.000 


215 


4+ 


2+ 


E2 


-0.332(1) 


-0.107(1) 


0.000 



* For this transition all \Jx = 1. 



C. The parity non-conserving E2/M2 multipole mixing ratio of the 501 keV tran- 
sition in i^°Hf™. 

The asymmetry, A, as defined by Eq. HI was evaluated for each field reversal (except 
5 and 7) for all transitions observed (the ^^Co transitions were not observed when the 
lead absorbers were present) with the results given in Table HI They show that there is a 
clearly established non-zero result of between 0.9% and 1.5% for the 501 keV transition in 
every reversal in which the nuclear sample was polarized. All other transitions show zero 
asymmetry within statistical error, and there is no evidence of asymmetry outside error in 
these and the 501 keV transition in the 'dummy' reversal 8. The larger experimental error 
for the 137 keV transition is caused by its lower intensity and poorer peak-to-background 
ratio than the other transitions (see Fig. [3]). 

Accepting the evidence that the 332 keV and 444 keV transitions show no effect, data 
for the 5th and 7th reversals were analysed using the ratio of the 501 peak counts to the 
sum of the 332 and 444 keV peak counts. It is shown in Appendix A that this ratio can be 
used to calculate the asymmetry of the 501 keV transition on the assumption that the other 
transitions have zero asymmetry. 

The asymmetries measured for the 501 keV transition are given in Table HIT] and plotted 
versus 1/T in Fig. [8] where they may be compared with theoretical calculations of the 
asymmetry obtained, using the same hyperfine interaction and fraction in good sites as 
mentioned above, for selected values of the E2/M2 mixing ratio e. The weighted average 
value of this mixing ratio is e = - 0.0324(16) (Table HTB . 
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D. Error analysis. 

The uncertainties in this experiment are of three types; statistical, theoretical and geo- 
metrical. The result given at the end of the previous section shows only the statistical error 
associated with the spectral data. This is simple to estimate as windows were set on the 
photopeaks, with no peak fitting, and background was taken from adjacent windows above 
and below the well-resolved peaks.. The other two sources of uncertainty are discussed in 
this section, leading to the final result. A second possible source of error in the 501 keV 
photopeak count is pile-up of pulses from detection of coincidences, either true or accidental, 
between 444 keV and 57 keV quanta. This effect has been estimated by comparing the weak 
pile-up peaks observed for other energy pairs (215 keV + 332 keV etc) with the product of 
their individual photopeak counts. This established that the contribution of pile-up events 
in the 501 photopeak was 0.4(1)% of the true single 501 quantum counts when the lead 
absorbers were absent, and orders of magnitude less when they were in place. Since, on field 
reversal, the change of total count rate is of order 1%, hence a maximum 1% change of the 
pile-up rate itself, pile-up has negligible influence on the measured asymmetry of the 501 
kev transition. 

The theoretical uncertainty derives from the adopted values of three parameters; the 
hyperfine interaction strength taken to be 25.9(13) mK, the fraction in good sites / measured 
to be 0.805(10) and the normal E3/M2 mixing ratio 6 in the 501 keV transition, taken to be 
5.3(3). To estimate the consequence of the uncertainties in these parameters on the value 
of the E2/M2 mixing ratio calculations of the asymmetry ^ as a function of e and inverse 
temperature, as in Fig. [H]were made using extreme values and compared with the standard 
calculations in which the central values, 25.9 mK, 0.805 and 5.3 were taken. The results of 
these calculations showed that uncertainty in the interaction, the fraction / and the E3/M2 
mixing ratio produced changes in values of e deduced from the measured asymmetries by, 
respectively, ± 1.9%, ± 2.7% and ± 4.9%. Added in quadrature, the total uncertainty, from 
these causes, in the average value of e is ± 0.0017. 

The geometrical uncertainty concerns the accuracy in positioning of the detectors, e.g 
their angles to the orientation axis, and in estimation of their solid angle correction factors. 
Allowing for a range of ± 2° in angle and of ± 1% in the solid angle correction factors 
produced changes in extracted e of about 1%, but in opposite directions. These small and 
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TABLE III: Results of asymmetry ^(T) in the 501 keV transition for all reversals, in order of 
decreasing temperature, and the extracted E2/M2 mixing ratio e. 



Reversal No. 


T[mK] 


1/T[K"^1 


A[%] 


e\%] 


6 


57(7) 


17(2) 


-1.48(26) 


-3.8(7) 


5 


25 Qd^) 


40(2) 


-1.19(13) 


-3.0(3) 


4 


22.2(14) 


45(3) 


-1.41(23) 


-3.7(6) 


3 


20.8(16) 


48(3) 


-1.12(13) 


-3.0(5) 


2 


19.2(13) 


52(4) 


-1.24(19) 


-3.4(5) 


1 


15.6(2) 


64(1) 


-1.09(32) 


-3.1(9) 


7 


7.6(1) 


132(2) 


-0.93(13) 


-3.4(5) 



Average e = - 3.24(16)% 



cancelling corrections have been omitted from the final error calculation. The high degree 
of cancellation between the A = 1 term, which is negative and the A = 3, 5 terms which are 
positive, in the asymmetry calculation, renders the result less sensitive to the solid angle 
corrections than might be expected and also insensitive to temperature except at very low 
values of 1/T. 

The final result is therefore e = - 0.0324(16) (17) where the first uncertainty is based on 
measurement statistics and the second the result of uncertain input parameters. 

VI. DISCUSSION. 

The measurements were designed to eliminate several possible sources of systematic er- 
ror. As detailed in Table [U reversals were carried out from different initial field directions, 
counting rates were varied by the addition of absorbers and reversals were done with and 
without the pulser, the presence of which gave rise to some peak distortion in the gamma 
spectra. The excellent stability of the beam spot position on the iron foil is demonstrated 
by the steadiness of the 332 keV/444 keV L/R ratio (shown in Fig. H]) throughout the exper- 
iment and also by the close equality of the 'warm' ratios taken at the start and finish of the 
experiment. The temperature was varied as widely as possible to give considerable change 
in the degree of polarization of the hafnium nuclei, thus providing the broadest possible 
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FIG. 8: Measured asymmetry A of the 501 keV transition, as a function of inverse temperature, 
compared with calculations using Eq. IHfor a range of values of the E2/M2 mixing ratio e. 

evidence for the temperature variation of the asymmetry. The results show the required 
robust self-consistency under all conditions. The final result for the irregular E2 /M2 mixing 
ratio in the 501 keV transition, e = - 0.0324(16) (17), is in extremely close agreement with 
the long accepted best value -0.030(2) Q|. Thus the present work, done with up-to-date and 
very different technology from the previous measurements, endorses fully the older results 
and has extended them. The observed temperature variation of the asymmetry, taken to 
significantly higher degrees of nuclear polarization than was accessible previously, shows full 
agreement with the behaviour predicted by the earlier result for e. This is a very satisfac- 
tory outcome which, if it lacks excitement, preserves the status of the asymmetry of this 
transition as the best established demonstration of parity admixture in nuclear phenomena. 
Theory thus faces the problem of providing understanding of the numerical result. Although 
the present experiment was not in every way fully optimized, there is little incentive to go 
further experimentally until some theoretical stimulus arises. 

Unfortunately, to date the problems presented by the specific nuclear structure of the 
levels involved, notably the high degree of K-forbiddeness of the 8~ decay, have prevented 
effective theoretical calculation of the expected degree of parity non-conservation in this 
case. Thus the very factor which is in all likelihood responsible for the large magnitude 
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£ = 


-0.028 


e = 


-0.031 


£ = 


-0.034 


£ = 


-0.037 




-0.040 




-0.044 




of the effect, the strong hindrance of the normal decay matrix elements, is also the cause 
of difficulty in giving a good theoretical description of the observed large parity violating 
effect. 
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VIII. APPENDIX A 

For reversals 5 and 7 normalization of the 501 keV transition to the pulser was not 
available and the ratio of the intensity of this transition to the sum of the combined 444 keV 
and 332 keV transitions was used as a measure of the 0°-180° asymmetry. Here it is shown 
that this is equivalent to the pulser normalized asymmetry provided the other transitions 
contain negligible parity violating terms in their angular distributions. 

For each transition and for a single detector at angle 9 to the axis of orientation, the 
recorded gamma count in any file, Ni can be written as the sum of product terms as in 
Eq. [H but it is convenient to separate the even, parity conserving terms from the odd, 
parity non-conserving terms, thus 

SfAW{e,El;)=Sc,{l + e, + o-:) (7) 

where e(o)i = BxUxAxQ xPx{cos 9) , S is the instantaneous source strength, Cj is the 

keven(odd) 

product of the fractional gamma emission per decay fi and the detector efficiency Si for that 
energy. 

The ratio R of counts in the 501 keV peak N501 to the sum of counts in the 444 keV peak, 
N444, and the 332 keV peak, N332, measured in the detector in the direction of the applied 
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field 9 = 0°(180°) pre(post) reversal, can be written: 

-Rpre(post)(^) ' 



^^"332 + ^444/ pre{post) 

5c50l(l + 6501(6*) + 0501(6')) 



(8) 



5C332(1 + e332(^) + 0332(^)) + ^0444(1 + e444(^) + 0444(^)) 

The difference between the ratios R measured in the initially 0° detector pre and post 
reversal, divided by their sum, after cancellation of the source strength, is given by 

^ _ Rpvc — -Rpost 

Rprc ~l~ -Rpost 

where the change of detector position from 6 = 0° (pre-reversal) to 9 = 180° (post-reversal) 
reverses the sign of the odd terms but does not change the even terms. 

After straightforward but somewhat lengthy manipulation this ratio becomes 

^ _ O50l[(c332(l + 6332) + £444(1 + 6444)] - £3320332(1 + 6501) - 04440444(1 + 6501) ^^Q^ 
£332(1 + e332)(l + 6501) + £444(1 + e444)(l + £501) - O501 (£3320332 + £4440444) 

When the product of odd terms (that is of parity non- conserving terms in the 501 keV 
transition and in either of the other transitions) in the denominator is neglected as small, 
this further simplifies to 

^ _ 0501 £3320332 + £4440444 ^-j^-j^^ 

(1 + 6501) £332(1 + 6332) + £444(1 + 6444) 

which can be seen to be half the asymmetry of the 501 keV transition minus half the weighted 
sum of the the asymmetries of the 332 keV and 444 keV transitions. 

If it is assumed that the other asymmetries are indeed zero, then the ratio TZ is simply 
half the asymmetry of the 501 keV transition. This assumption is adopted in the analysis of 
reversals 5 and 7, both because the experiment provides evidence for the other asymmetries 
being small and because of the theoretical difference between the strong intra-band nature 
of the 332 keV and 444 keV E2 transitions and the highly K-forbidden 501 keV M2/E3 
transition. 

In reversals 5 and 7 the two axial detectors yield independent values of TZ, TZ{0°) and 
7^(180°), of opposite signs, since one detector is initially at 0° and the other at 180°. The 
combination 7^(0°) - 7^(180°) (equals 2 72.(0°)) gives the full asymmetry of the 501 keV 
transition. This is given in Tables [T] and IIIII 
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